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We analyze the valley composition of one electron bound to a shallow donor close to a Si/barrier 
interface as a function of an applied electric field. A full six- valley effective mass model Hamiltonian 
is adopted. For low fields, the electron ground state is essentially confined at the donor. At high 
fields the ground state is such that the electron is drawn to the interface, leaving the donor practically 
ionized. Valley splitting at the interface occurs due to the valley-orbit coupling, V^ = \Vj \e t6 . At 
intermediate electric fields, close to a characteristic shuttling field, the electron states may constitute 
hybridized states with valley compositions different from the donor and the interface ground states. 
The full spectrum of energy levels shows crossings and anti-crossings as the field varies. The degree of 
level repulsion, thus the width of the anti-crossing gap, depends on the relative valley compositions, 
which vary with \Vj Q \j and the interface- donor distance. We focus on the valley configurations of 
the states involved in the donor- interface tunneling process, given by the anti-crossing of the three 
lowest eigenstates. A sequence of two anti-crossings takes place and the complex phase affects 
the symmetries of the eigenstates and level anti-crossing gaps. We discuss the implications of our 
results on the practical manipulation of donor electrons in Si nanostructures. 
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I. INTRODUCTION 

The search for a functional quantum computer (QC) 
started in the mid-nineties, and by the year 2000 many 
systems had been considered as candidates for its phys- 
ical implementation.^ Among them, the 1998 proposal 
for a Si-based QC by Kane 2 raised special interest due to 
objective and relevant factors favoring Si, such as the ac- 
cumulated know-how in processing Si for advanced device 
applications, the relatively long spin coherence times and 
the possibility of isotopic purification processing, further 
increasing coherence times. 3 

On the other hand, the conduction electrons in Si are 
not in a well defined single Bloch state. Instead, the 
Si conduction band is six-fold degenerate, with minima 
(valleys) along the x, y and z crystallographic directions. 
This imposes limitations to the spin manipulation and 
coherence)^ 

It was recently proposed to encode quantum infor- 
mation directly into the valley degree of freedom, con- 
verting the spurious valley Hilbert subspace into a use- 
ful ingredient for a QC. 5 Naturally, this raises funda- 
mental questions, such as how to promote controlled 
inter- valley transitions, to what extent valley degener- 
acy can be lifted, and how sensitive such operations are 
to fabrication-related parameters. The valley degree of 
freedom also affects transport properties in Si nanostruc- 
tures: valley degeneracy has been recently shown to pro- 
duce a valley Kondo effect in a singly doped Si fin field 
effect transistor. 6 

We study here the valley degree of freedom for one elec- 
tron bound to a donor — more specifically substitutional 
P in Si — tunnel-coupled to a (001) Si/Si02 interface at 



a distance d from the donor. The barrier material is 
taken to be Si02 for definiteness, but it could in princi- 
ple be any high quality interface, such as Si/SiGe. The 
evolution of the inequivalent valley contributions is ob- 
tained by mapping the low-lying manifold (following in 
more detail the three lowest energy states) as an elec- 
tric field pulls the electron away from its "hydrogenic" 
configuration at the donor site towards the state at the 
interface. 

A ground state electron confined in the direction per- 
pendicular to the interface in the triangular potential 
formed by the barrier and the electric field [see Fig. HJ a)], 
still remains bound to the donor core potential, leading 
to localization in the in-plane direction. 78 For low fields, 
the electron in the ground state is essentially confined at 
the donor, where the lowest energy manifold is split by 
the tetrahedral crystal field environment into states with 
distinct contributions from the six valleys. In particu- 
lar, the ground state at the donor is a non-degenerate 
symmetric combination of the six valleys. 

At high enough fields, the donor is ionized, and its 
electron is shuttled to the interface. At the interface, 
the valley levels split into a four-fold degenerate excited 
state, consisting of ±x and ±2/ Bloch states, and a lower 
manifold spanned by the ±z valleys. The two lowest lev- 
els are only slightly separated due to the abrupt interface 
breaking the z reflection symmetry (generally less than 
1 meV separatiorP^S) . 

These relatively simple and well understood valley 
compositions at interfaces and isolated donors could lead 
to a variety of compositions at intermediate fields, as 
illustrated by the main panel in Fig. [I] The lines are 
symmetry-allowed paths connecting energy levels from 



the low-field (right of the panel) to the high-field (left) 
regime: one may anticipate a rich variety of behaviors 
and formation of hybrid donor-interface states. Each 
level crossing or anti-crossing as the field varies may 
change the valley compositions of the involved states. 

The Stark shifted spectrum of P donors in bulk Si (no 
barrier material) is discussed in Ref. [15]. The effective 
mass approach is adopted and the envelope functions are 
expressed as a combination of atomic-like orbitals. The 
same group considered later the Stark effect for P donors 
at the center of Si nanospheres embedded in a barrier 
material.^ The quasi-spherical symmetry of the problem 
without field allows the assignement of hydrogenic quan- 
tum numbers to the interface states. Tight-binding cal- 
culations for systems more directly related to the present 
geometry were reported in Ref. [17]. In these references 
elaborate numerical procedures are adopted leading to 
accurate results. Here we get good qualitative agree- 
ment for the spectrum in comparison with these pre- 
vious works, and we explore a complementary aspect - 
namely the valley-composition evolution of the low-lying 
states under an increasing external field. The same ge- 
ometry has been previously studied in a two- valley model 
{(000001), (000010)} where the valley-composition anal- 
ysis is not accessible.^ 

We discuss here the evolution of the valley quantum 
number as a function of the electric field for field values in 
the range where the hybridization between donor and in- 
terface states takes place. We show that the phase of the 
interface valley-orbit coupling plays an important role in 
defining the gap amplitudes, affecting manipulation ca- 
pabilities for the donor electrons. Our model sheds light 
to the qualitative features of the level diagrams shown 
in Refs. |15|17| , and permits to predict the expected dia- 
grams for various geometries of the donor/barrier prob- 
lem. Our aim is to get a clear description of the changes 
in symmetry and valley composition of the states in- 
volved in the donor-interface electron shuttling problem. 
Our results may be useful for donor electron valley ma- 
nipulation via an applied field. 

This paper is organized as follows. In Sec. [IT] the sys- 
tem is described and a 6-valley effective mass model- 
Hamiltonian for the donor electron under an applied elec- 
tric field is explained. In Sec. |III[ our results for the 
electric field dependence of the low-lying spectrum and 
of the valley composition of the three lowest electronic 



eigenstates are presented. We conclude in Sec. IV with 
a discussion of the implications of these results on the 
practical manipulation of donor electrons in Si nanos- 
tructures. 



II. MODEL 

We consider a single electron bound to a substitutional 
P donor at z = near a Si (001) / Si02 interface at 
z = —d and under an applied uniform electric field per- 
pendicular to the interface, F = Fz pointing from the 
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FIG. 1: (Color online) (a) Double well potential in the z- 
direction formed by the Coulombic donor potential plus the 
triangular interface/electric field potential, d is the distance 
between the donor and interface. (Main panel) Symmetry of 
levels at the donor and at the interface. Every level is de- 
scribed by six coefficients corresponding to the six valleys of 
Si conduction band: (x,—x,y,—y,z,—z). This defines the 
valley composition of each state. At the interface, the mass 
anisotropy breaks the valleys degeneracy in a doublet (z,—z) 
and a quadruplet (x,—x,y,—y). The doublet degeneracy is 
lifted due to the valley orbit coupling (Vj = \Vj \e l ) aris- 
ing in a sharp (001) interface, 1 9 1 12 ' 14 1 as shown in (b). Cf and 
Cj z are defined in Eq. (I6J), and Cj z — —Cj z . At an iso- 
lated donor, the valley orbit coupling leaves a no n- degenerate 
ground state with A\ symmetry, well separated (splitting ^12 
meV) from the other five levels ! 20 " 22 ' The binding energies 
given on the right of each level are experimental values for 
bulk P donors" The lines join the valley compositions at 
donor and interface that are connected by symmetry. 



barrier to Si, thus pushing the electron away from the 
donor and towards the interface (see Fig. W[ a)). In effec- 
tive "atomic" units for Si, Ry* = m_ L e 4 /2fi 2 e| i = 19.98 
meV and a* = h 2 esi/m±e 2 = 3.157 nm, the Hamiltonian 
is written a d 8 * 19 1 
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d 2 



7- 
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dx 2 dy 2 dz 2 



keFz + H v 



(1) 



with 7 = raj_/ra|| as the ratio between the trans- 
verse (m± = 0.191m) and longitudinal (ray = 
0.916m) effective masses, es% = 11.4 , k = 3.89 • 
10 _7 e| i (m/m_L )cm/kV, and the electric field F is given 
in kV/cm. In Eq. (M), the kinetic energy is { — -§^2 — 
Jf-2 — 7^2)5 the next two terms are the donor Coulomb 
potential and the electric-field linear potential, respec- 
tively. The last term describes the valley-orbit effects, 
namely the coupling between different valleys due to the 
singular nature of both the donor (D) and the interface 
(I) potentials, as described below. 

As shown in Figflla), the system can be modeled by 



the combination of two potential wells: one that binds the 
electron to the donor (at low fields) and another binding 
the electron at the interface (at high fields). 8 Note that 
the interface potential includes not only the electric field 
perpendicular to the interface but also the Coulomb at- 
traction to the donor at a distance d which confines the 
electron in the xy-plane even when it is at the interface. 7 
We simplify the calculation by initially computing the 
variational ground state wave-functions at the donor and 
at the interface, which define the envelopes in the basis 
for the full six valley problem. 

The conduction band of Si has six degenerate minima 
(valleys) in the (100) directions at a distance fco = 2tt^^ 
from the V point, where asi = 5. 4 A is the lattice param- 
eter of Si. In the simplest effective mass approximation, 
only the Bloch functions at the positions of the conduc- 
tion band minima are considered, and the ground state 
of the electron at the donor is written 20 



^d = 



J2 C£F£(r)^(r), 



(2) 



/X = ±£C,±2/,±2 



r D 



(r) are envelope functions and 
are the six Bloch eigenstates at the 



where Fj^ (r) 

conduction band minima. We take the variational donor 
envelope functions Fj^{r), where v = |/i|, centered at 
r = 0, following the form introduced in Ref. [20] 
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These are normalized hydrogenic Is envelopes, distorted 
due to the Si conduction band effective mass anisotropy, 
and {AT]}} are normalization factors. The effective Bohr 
radii a and b are variational parameters chosen to mini- 
mize the ground state energy. For the distances used here 
(d > 2a* ~ 6 nm), a and b coincide with Khon and Lut- 
tinger's variational parameters for a single impurity in 
the bulk 20 (d —> oo), namely, a = 2.365 nm and b = 1.36 
nmP 

The six-fold degeneracy of the ground state is lifted at 
a substitutional impurity because the translational sym- 
metry of the host crystal is broken, leading to interval- 
ley scattering effects known as the valley-orbit interac- 
tion) 2 ^ 2 ^ This effect can be accounted for phenomenolog- 
ically, introducing a coupling between valleys in perpen- 
dicular directions (e.g. x, z) — A c and in parallel direc- 
tions (e.g. z,—z) — A c (l + S c ) 25 This splits the unper- 
turbed six-fold-degenerate donor electron ground state 
into a singlet (A\ symmetry), a triplet (T2 symmetry) 
and a doublet (E symmetry), see Fig. [Tl For P in Si, the 
relative splittings between the different symmetry levels 
are reproduced taking A c = 2.16 meV and 5 = — 0.3™ 

The degeneracy is also lifted near the (001) 
interface! 26 ^ 27 ! First, due to the mass anisotropy, the z 
and — z perpendicular valleys are lower in energy than 



the ones parallel to the interface. The two lowest energy 
states are combinations of the z and — z valleys, whose 
double degeneracy is lifted due to valley-orbit coupling 
Vy Q at an abrupt interface. In general, Vj is a com- 
plex quantity with an absolute value proportional to the 
applied electric fielcP^ |V^| = AF, and dependent on 
the barrier height and abruptness.^ The prefactor A has 
been estimated by several authors.^ We use initially the 
largest, A = 1.36 A, as suggested in Ref. [11 . For in- 
stance, for F = 50 kV/cm, we have |V^| = 0.68 meV. 
We also consider the smaller value of A = 0.215A, esti- 
mated by Sham and Nakayama, 9 and discuss the quali- 
tative changes that occur in the spectrum. The complex 
phase of this valley-orbit coupling is also dependent on 
the interface quality and has been estimated to be ~ — 1 
for an abrupt Si/Si02 interface.^ H vo in Eq. (fTl) takes 
into account all these valley-orbit interactions, both at 
the donor (D) and at the interface (I). Its form (in a 
basis set defined next) is described in the Appendix. 

Following Eq. (pi), we take the lowest states at the in- 
terface as 



#1 



/j,=±x,±y,±z 



C?F?(r)Mr)- 
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The envelope functions -Fj (r), with v 
the variational forrrP 



|/i|, are taken in 



Ff = Nf(z + d) 2 e- ax y (z+d)/2 e-^ lX +(32y )/2 , (7) 

F V j = 7V|( z + ^)2 e -^(^)/2 e -(/3 2 . 2 +/3 iy 2 )/2 ? (g) 

(9) 



Ff = 7Vf(z + d) 2 e-^ +d )/ 2 e-^/ 2 , 



where the penetration into the barrier is considered to 
be negligibly small. Here, p 2 = x 2 + y 2 , and a, a xyi ft, 
Pi and /?2 are variational parameters: 1/a and l/a xy are 
related to the width of the wave functions along the z- 
direction, and depend on the value of the applied electric 
field; /?, Pi and P% correspond to the confinement in the 
xy plane, which is controlled by the attractive potential 
of the donor, hence the "/?" parameters depend on the 
distance dP 

In Eqs. Q and (J6J) we do not include explicitly the 
pinning point for the plane- wave part of the Bloch func- 
tions at the the donor and interface potentials. Effects of 
the interference induced by the different pinning points 
are discussed in Ref. [18]. We solve for the lowest states 
of the full potential (donor, electric field and interface) 
by obtaining the spectrum of H in the combined basis 
of the lowest D and I envelopes as determined variation- 
ally, each multiplied by the respective Bloch functions. 
This defines the Hilbert space for our model calculation. 
States with different Bloch indices \i 7^ v are not cou- 
pled unless there is a non-zero contribution from H vo (see 
Appendix A). Within the Hilbert space defined here, the 
Hamiltonian is represented by a 12 x 12 matrix, written 
formally as four 6x6 blocks 



H 



Hdd Hid 
Hdi Hjj 



(10) 



The equation giving the spectrum takes into account the 
non-orthogonality of our basis, i.e., we solve for H^i — 
EiS^i, where S is the 12x12 overlap matrix so that 
S^d = S^ = <W and S% = S% = 5, v (F£ \ F?>. 
Here, (F£ \ Fj) is the overlap between the interface and 
the donor envelope functions, an exponentially decreas- 
ing function of d. A general state is written as: 



*= E E 

H=±x,±y,±z L=I,D 



C£F£(r)0 M (r). 



(11) 



The complex coefficients {C£} give the contributions to 
the state \P of each valley /i at the interface or donor 
(L = I,D). Also, they define the symmetry of the state, 
and are referred to here as "valley coefficients". The 
normalization condition reads 



X) ci*sr v c\ 



Era 2 



cross terms = 1, 



(12) 
where the cross terms are zero for L = V and/or for 
\i 7^ v. For d = 4a* the cross terms due to the finite 
overlap between donor and interface envelopes are less 
than 5% in the electric field range of interest. The overlap 
Sid is negligible for d = 5a* and larger. 

In Fig. [I] we show the valley coefficients of the states at 
an isolated donor (just the {C^} are given since Cj = 
for all ji) and at an interface under a perpendicular elec- 
tric field (just the {Cj } are given since C^ = for all fi). 
The valley coefficients are perturbed when the impurity 
is located at a distance d from the interface, and hybrid 
states may be formed P^ where both D and / coefficients 
contribute. The different initial (for low F) and final (for 
strong enough F) states are symmetry-compatible when 
their valley compositions are not orthogonal [see Figjl] 
(main panel)]. 



III. RESULTS 



Previous theoretical and experimental studies^- 17 * 29 * 30 * 
identified and analyzed a characteristic field at which 
the electron ground state crosses over from being bound 
mainly around the donor (donor-like) into being bound 
mostly near the interface (inter face- like). This field, F c h 
(see Fig. [2|, decreases with the distance d from donor 
to interface, while the tunneling time increases exponen- 
tially with d. These properties are consistent with our 6- 
valley model results. We explore here the valley degrees 
of freedom, mainly close to the characteristic field, and 
analyze how the applied field changes the valley compo- 
sition of the main electronic states involved in the donor- 
interface shuttling. 

The electric-field dependence of the complete spectrum 
of the Hamiltonian in Eq. (10), for d = 4a*, is shown in 
Fig. [2] The overall level structure here is similar to the 
one presented in Refs. |15|17j . At small fields, F < 32 
kV/cm, the six lowest eigenstates correspond to donor- 
bound states. The degeneracies of the T2 triplet and 
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FIG. 2: (Color online) Full spectrum of eigenvalues for d = 
4a* in a wide range of electric fields. For small values of 
the electric field, the lowest six eigenvalues correspond to the 
donor states and the highest six to interface states. We use 



here V' = IK, 



Z7T/3 



with IKoI = XF and A = 1.36 A. 



the E doublet are lifted (with a very small splitting) 
due to the perturbation produced by the interface. At 
large fields F > 70 kV/cm, the six lowest eigenstates be- 
come interface-bound states. The two lowest eigenstates, 
which are combinations of the z and — z valleys, are split 
by 2|V^|. Different eigenvalues cross over each other as F 
changes. The levels cross or anticross depending on their 
relative symmetry. The general scheme shown in Fig. [2] 
is qualitatively similar for different values of the distance 
d. The size of the gap at the anti-crossings is related to 
the shuttling time and decreases as d increases. 8 

We follow now on the evolution of the valley contri- 
butions of the three lowest eigenstates, which, via mu- 
tual energies crossing over, "become" or contribute to 
the ground state for some range of field values (see Fig. 
[3J) . Because of time-reversal symmetry, valleys \i and — \i 
contribute equally to any given eigenstate, so we may 
quantify the valley contributions by valley populations 
of each direction v — x, y, z, defined as 



H=±v L=I,D 



(13) 



Due to the cylindrical symmetry of the system, P z is in 
general different from P x and P y , while the two latter 
are equivalent. This means that states with exchanged 
coefficients CjfL 1 d ^ ^l^i d &i ye ^ ne same expectation 
value for the energy, since H is invariant under ±x <-} ±y. 
So P x and P y are presumably equal. It may occur for 
a particular state that the weights P x and P y obtained 
numerically differ: in this case a degenerate state is al- 



ways found with P' x 



P y and P y 



P x , as expected. 



Therefore, differences in P x and P y are not physically 
meaningful and we present our results in terms of P z and 



(P x + Py)/2 



£ £ m 2 . (14) 

=±x,±y L=I,D 



In this definition we do not take the cross terms from 



Eq 



• Pfi 

allowir 



into account, so normalization gives 2P X 



1, allowing P z and P xy to be directly compared to each 
other, giving the relative weight of the z and the aver- 
age x and y populations. The lowest interface state at 
large F only involves z and — z valleys, therefore P xy = 
and P z = 1. On the other hand, the lowest donor state 
involves a symmetric combination of all valleys leading 
to P xy = P z . Hybrid states correspond to intermediate 
values of P xy and P z : < P xy < 1/3 and 1/3 < P z < 1. 
In Fig. p] (upper frames) the spectrum is presented for 
a reduced range of energy and fields around the charac- 
teristic field for (a) d = 4a* (F ch ~ 46 kV/cm) and (b) 
d = 5a* (F c h ~ 35 kV/cm). The three rows of frames 
below give the corresponding valley populations for the 
2 nd excited, 1 st excited, and ground (GS) states, respec- 



tively. Here we take a complex V^ 



IK 



t/3 



From 



Fig. [3] one can clearly observe that the anti-crossing at 
F c h in fact involves two anti-crossings: one between GS 
and 1 st excited, and another one between 1 st and 2 nd ex- 
cited states. Well below F c h, the GS is donor-like with 
Pxy — Pz- Above the two anti-crossings the GS and 1 st 
excited are interface-like states, only involving the z and 
— z valleys, thus P xy = 0. 

We note that the ground state and the 2 nd excited 
state swap their valley compositions for fields below and 
above the crossover region. In fact, comparison of differ- 
ent frames in Fig. [3] show that the low-field D-like com- 
position in (g) "moves" to the high- field behavior in (c), 
while the low-field I-like composition in (c) is found in 
the high- field behavior in (g). Along the crossover re- 
gion the compositions change smoothly or abruptly (ac- 
cording to d) among the limiting behaviors. The same 
applies to frames (d) and (h) at a more distant donor po- 
sition from the interface. The intermediate state, shown 
in (e) and (f ) , preserves the compositions at low and high 
fields, while it is clear that this state hybridizes with both 
ground and 2 nd excited states along the crossover region. 
In summary, the net effect of the field on the three lowest 
eigenfunctions far from F c h is to cross the ground and 2 nd 
excited states, while the intermediate I s * excited state is 
not affected overall, although all three states mix at the 
crossover range. An extreme example of the I s * excited 
state not being affected by the crossover is discussed be- 
low, in the context of FigJ5] 

Comparison between (a) and (b) in Fig. [3] illustrates 
the well known effect of increasing d, already discussed 
in previous publications pOD namely reducing anti-crossing 
gaps and sharpening transition lines. Between the two 
anti-crossings, the 1 st excited state is donor-like for d = 
5a* in Fig. |3|f), however for d = 4a* in (e) the 1 st excited 
state is an interface-donor hybrid due to the larger donor- 
interface overlap. We note that all the gaps obtained for 
d = 5a* are extremely small compared to the relevant 
energy scales here, so for most practical purposes the 
behavior would be equivalent to level crossing. 

The phase on Vj Q = | V^|e 26> , not as extensively con- 
sidered so far, is analyzed in Fig. [4] for d = 4a*. The 
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FIG. 3: (Color online) Evolution of the valley population for 
the three lowest eigenvalues around the characteristic field 
(see Fig. [2]) for two different distances (a) d — 4a* and (b) d — 
5a*. We use Vj Q = |V£|e i7r/3 with |v£J = XF and A = 1.36 
A. The top panels reproduce the eigenvalues involved in the 
lowest energy anti-crossing, which is in fact a sequence of two 
anti-crossings. All other panels show the valley population 
in the second excited state, the first excited state, and the 
ground state (GS) in different lines. The red (solid) curves 
correspond to the weight of the ±2 valleys or longitudinal 
weight (at donor and interface), and the green (dashed) curves 
are the weight of the ±x and ±y valleys or transversal weight 
(at donor and interface). Labels D or I refer to donor- like or 
interface- like states in terms of real space location. Here D 
is a combination of the 6 valleys and I involves the z and -z 
valleys. 



complex phase of the valley-orbit coupling at the in- 
terface affects the symmetries of the eigenstates leading 
to different gaps at the two anti-crossings around F^. 
Figs. J4^a) and (g) show the limiting cases of = and 
= 7T, which correspond to a real Vj Q , lead to a zero-gap 
(crossing) involving the symmetric donor-like eigenstate 
and the antisymmetric interface state, with all C£ = 0, 
except Cf = —Cj z = 1/a/2. For a general 0, the two 
lowest interface states always have a symmetric compo- 
nent which "repel" the symmetric donor-like level. The 
size of the gap at the anti-crossing increases as the weight 
of the symmetric part of the interface- like state becomes 
larger. The two gaps become equal for 6 = 7r/2 which 
corresponds to a purely imaginary Vj Q . In summary, for a 
fixed I Vy Q I and calling gz, and qr the gaps to the left and 
to the right in energy, we get g^ > gn for < < tt/2 
and qr > gL for tt/2 < 6 < tt. In particular g R = 
(g L = 0) for = (tt) and g L = g R for = tt/2. The 
largest gap observed for d = 4a* is ~ 1 meV, of the same 
order than the interface valley-orbit splitting considered. 
In contrast, for d = 5a* the gap is ~ 0.035 meV, almost 
two orders of magnitude smaller. 

The results presented so far correspond to a relatively 
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FIG. 4: (Color online) Three lowest eigenvalues around the 
anti-crossings region, corresponding to d = 4a* for different 
values of the phase of the valley-orbit coupling at the in- 
terface (Ko = \Vj o \e i0 ). |Ko| = AF and A = 1.36 A. The 
extra panel on the right top corner shows the values of the 
two anti-crossing gaps versus 0. The dashed line represents 
the gap for the first anti-crossing between the GS and the 
1 st excited state while the solid line is the gap between the 
1 st and 2 nd excited states. Note that for 6 = and 7r, Vj Q 
is a real quantity and one of the anti-crossings has zero gap 
(namely, it is actually a 2- level crossing). The two gaps are 
equal for = tt/2 which corresponds to a pure imaginary Vj Q . 
Our results are obviously invariant for <H> —0. 



large value of the valley-orbit coupling at the interface, 
with A = 1.36 A as estimated in Ref. [11 J. For this case 
and d = 4a*, the gap at anti-crossing is slightly smaller 
than the value of the valley-orbit splitting (2|V^|). Pre- 
vious calculations by Sham and Nakayama 9 lead to a 
smaller A = 0.215A. In Fig. [5J the three lowest eigenval- 
ues close to the characteristic electric field for d = 4a* 
for the Vj Q as calculated by Sham and Nakayama are 
shown. Here, the gap at anti-crossing (at F = F c h) is 
about three times larger than the value of the valley- 
orbit splitting at the interface for d = 4a*. Due to the 
relatively smaller value of the valley-orbit splitting, the 
interface states are much closer: The two anti-crossings 
at F c h seem to merge into a single one and the 1 st ex- 
cited state is always an interface state. This is in con- 
trast with results for the larger value of V^ Q where the 
1 st excited state is hybridized or donor-like between the 
two anti-crossings [see FigJ3le) and (f)]. Another distinct 
feature of the small \Vj Q \ limit, illustrated in Fig. ^1 [com- 
pare with Fig. |4jc)], is the way the three levels separate 
around F c h, with level repulsion among the outer ones, 
and no deviation of the middle state from the linear path: 
it does not couple to the others. Fig. [5] is very similar 
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FIG. 5: (Color online) Three lowest eigenvalues for d = 4a* 
close to the characteristic field with the valley orbit coupling 
at the interface |V^,| — ^F with A = 0.215 A as calculated 
by Sham and Nakayama— We use = tt/3 as the phase of Vj Q 
[compare with Fig.CTc)]. 



to the inset of Fig. 2 in Ref. [15], where the model does 
not include a barrier, thus corresponding to the Vj = 
limit. We may infer from these results that different ion- 
ization regimes for doped Si may arise as a function of 
d (inversely related to the gap at anticrossing) and the 
valley-orbit coupling at interface (related to the electric 
field, to the height, and to the quality of the interface 
barrier) . 



IV. DISCUSSION 

We have analyzed in full depth the valley contribu- 
tions to the three lowest energy levels of a single elec- 
tron in the donor-in-Si-near-a-barrier system under an 
applied external field, focusing in the field range around 
the crossover between donor-like and interface-like char- 
acter of the ground state. Strong hybridization occurs in 
the vicinity of the crossover, which we identify as a se- 
quence of two anti-crossings (see Fig.[3|. A result to keep 
in mind in practical applications is the strong dependence 
of the levels crossings and anti-crossings on the phase of 
the valley-orbit coupling at the interface, Vj Q = l^ol 6 * 61 - 
This is to be expected since the phases affect directly the 
symmetry of the states. A trivial example is the case 
9 = 7T (0), where the GS is a symmetric (anti-symmetric) 
combination of the z and — z valleys. It is not a straight- 
forward task to predict or control the phase #, which 
should vary with the barrier materialp' interface rough- 
ness due to steps, inter diffusion, etc, and other sample 
properties. 

In practice, the results on the double anti-crossings 
shown in Fig. [1] may play an important role in applica- 
tions involving the donor /barrier system. For instance, 
in the limit of — >• 0, the ground state is well separated 
from the excited states, suggesting the possibility of adi- 
abatically shuttling the electron from the donor site to 
the barrier interface. This could be a suitable scenario 
for a spin qubit for which one needs the valley splitting 
to be larger than the Zeeman splitting. An intrinsically 
large |V^,|, like the one considered in Figs.[2j|3j and|4j is 
also needed, as well as a small d (< 4a*) to guarantee a 
sufficiently large gap at anticrossing. 



On the other hand, following the adiabatic theorem, 
the donor- interface shuttling would be prohibitively slow 
for 9 — >• 7r due to the crossing between the ground state 
and the first excited level. Another phase-related effect 
occurs due to the variable pinning point in the Bloch 
functions <^(r) |18|31j . which produces a periodic de- 
pendence of the valley-orbit splitting on d leading to the 
closing of the gap at anticrossings at some particular val- 
ues of d [H] . 

If instead of the spin we want to use the valley degree 
of freedom to define qubits, a controllable valley-orbit 
coupling is required.^ This is achieved if the energy sepa- 
ration of two states with the same envelope function but 
different valley compositions (therefore, different oscilla- 
tions in the atomic scale) varies significantly with the 
external applied fields. This situation can be attained 
above F c h when the two lowest eigenstates are mainly 
interface states involving different combinations of the 
z and — z valleys. A strong dependence of the splitting 
with electric field can be found close to F c h in some cases. 
For instance, in Fig. |4^g) (corresponding to 6 = tt) the 
closing of the gap between the two lowest eigenvalues pro- 
duces a fast decrease of the splitting as the electric field is 
lowered from ~ 48 to ~ 46 kV/cm. Also, for smaller val- 
ues of Vj , as exemplified in Fig. [5j the inter face- donor 
hybridization leads to a relatively large gap compared 
to Vy Q and consequently the splitting increases fast from 
~ 48 to ~ 46 kV/cm. However, the variation of the level 
splitting on the electric field is stronger when the donor- 
interface hybridization is large. This hybridization mixes 
valleys in the different directions, potentially producing 
decoherence in the valley sector. 

On the other extreme, = 0, Fig. |4^a) reveals that 
the upper valley composition crosses the second excited 
state: As a consequence, the valley information would be 
lost. 

In summary, we show that both the modulus and the 
phase of Vj Q affect the quantum behavior of donor elec- 
trons in Si near an interface. The value of is hard 
to predict, and its calculation probably requires knowl- 



edge of the atomistic distribution at the interface and 
a complete electronic structure description of the com- 
bined Si/structured interface/barrier system. It is not 
clear whether a direct experimental measure is possible. 
We propose that measurement of the separate gaps may 
provide an estimate of 0, as shown in the inset of Fig. [4] 
Finally, the trends and discussions regarding Fig(3]should 
bring new and valuable insight towards controlled valley 
manipulations. 
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Appendix A: Valley-orbit term (H vo ) 

The last term in Eq. 0, with A± = -A c = -2.16 
meV, Ay = -A c (l + S c ) = -1.51 meV, and \Vj Q \ = XF 
with A = 1.36 A. This is a 12x12 matrix in the basis set 
of the six valleys at the donor and the six valleys at the 
interface. According to the notation in Eq. (10): 
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